Linear E model in the Gaussian wave functional 
approximation II: Analyticity of the S-matrix and the 

effective potential/action 

Issei Nakamura^ 1 '' and V. Dmitrasinovic^ 
(1) Research Center for Nuclear Physics, Osaka University 
Ibaraki, Osaka 567-0047, Japan; and 
(2) Vinca Institute of Nuclear Sciences 
P.O.Box 522, 11001 Beograd, Yugoslavia 

February 1, 2008 

Abstract 

We show an explicit connection between the solution to the equations of motion 
in the Gaussian functional approximation [3] and the minimum of the (Gaussian) 
effective potential/action of the linear E model, as well as with the N/D method in 
dispersion theory. The resulting equations contain analytic functions with branch 
cuts in the complex mass squared plane. Therefore the minimum of the effective 
action may lie in the complex mass squared plane. Many solutions to these equations 
can be found on the second, third, etc. Riemann sheets of the equation, though their 
physical interpretation is not clear. Our results and the established properties of 
the S-matrix in general, and of the N/D solutions in particular, guide us to the 
correct choice of the Riemann sheet. We count the number of states and find only 
one in each spin-parity and isospin channel with quantum numbers corresponding 
to the fields in the Lagrangian, i.e. to Castillejo-Dalitz-Dyson (CDD) poles. We 
examine the numerical solutions in both the strong and weak coupling regimes and 
calculate the Kallen-Lehmann spectral densities and then use them for physical 
interpretation. 



1 Introduction 

In the present paper we extend the study of a chirally invariant, Lorentz invari- 
ant, self-consistent mean-field, variational approximation, that goes by the name of 
Gaussian wave functional approximation [1, 2] to the linear sigma model, that was 
begun in Ref. [3]. We have shown in Ref. [3] how to ensure chiral symmetry in the 
Gaussian approximation method, a major improvement over previous treatments. 
A number of questions have remained open after that paper, however. In particular 
we have not addressed the connection between the Gaussian approximation to the 
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canonical equations of motion and the Gaussian effective potential (EP) method, 
that is rather popular in finite temperature applications [4] . There the meson masses 
are defined in terms of the curvature (second derivative) of the effective potential 
evaluated at the minimum. It used to be believed that this definition leads to a vio- 
lation of the Nambu-Goldstone (NG) theorem for the tt fields [5], even in the chiral 
limit. This misunderstanding was cleared up in Ref. [6]; moreover it was shown 
there that one ought not to minimize the effective potential, which is momentum 
independent; rather one must minimize the effective action (EA) which leads to 
momentum-dependent equations. The latter mass definition leads to an equation 
to be solved for m a . Yet there have been no attempts to solve this equation in 
the literature, save for Ref. [3]. Thus, most m CT values present in the literature 
are not acceptable. The equation for m a is a transcedental one, however, with 
infinitely many Riemann sheets and an apparently indeterminate number and prop- 
erties of solutions. Furthermore, the Gaussian method involves certain auxilliary 
objects, such as the solutions to the gap equations (that are often interpreted as 
meson masses) whose physical role is also unclear. Similarly, the canonical Gaussian 
approach involves two-body scattering Bethe-Salpeter (BS) equations that do not 
seem to appear in the EP approach. 

In this paper we shall answer the aforementioned questions and some others not 
mentioned above: For example we show that exactly the same equations for the me- 
son mass appear in the canonical Gaussian approximation and the Gaussian effective 
action approach: indeed these equations are the net result of the coupled Bethe- 
Salpeter and the gap equations. Of course, this fact does not make them any easier 
to solve, but it offers a useful perspective on the number and nature of the solutions. 
Inhomogeneous BS equations are scattering equations, and in this particular approx- 
imation they will be shown to be equivalent to N/D equations of dispersion theory 
[7, 8, 9, 10], that ensure manifest unitarity. Some properties of their solutions, such 
as analyticity, and the physical interpretation of the solutions, follow from unitarity 
and causality, and have been known since the early 1960's [7, 8, 9, 10]. Another well 
known property of N/D equations is the arbitrary number of their solutions: this is 
the Castillejo-Dalitz-Dyson [CDD] ambiguity [7, 8, 9, 10]. The Gaussian approxima- 
tion is more restrictive than the N/D approximation, however: all properties, such 
as the number of CDD poles and values of subtraction constants are determined by 
the gap and BS equations that are a part of the canonical Gaussian approximation. 

In Ref. [3] we have numerically solved the BS equation in the scalar channel on 
the real s axis and found multiple solutions for certain parameter values, and no 
solutions at all for others. Yet, in the weak-coupling limit the Gaussian solution 
is unique and smoothly connected to the perturbative one. We use the Kallen- 
Lehmann spectral representation to show that the "heavy" solutions do not have 
particle- like properties. In consequence of this, we show that there is only one 
solution in each spin-parity-isospin channel. 

Thus we have established an explicit connection between some previously sep- 
arate formalisms, such as the effective potential/action (a.k.a. Cornwall-Jackiw- 
Tomboulis (CJT)) method [11], the Hartree + Random Phase Approximation (RPA), 
and the N/D method in dispersion theory, as well as shed light on the particle con- 
tent in this approximation. 

This paper falls into five sections. After the Introduction, in Sect. II we outline 
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the Gaussian effective potential/action method and prove its equivalence with the 
Gaussian approximation to the canonical equations of motion. In Sect. Ill we 
demonstrate the latter's connection with the N/D equations of S-matrix theory. In 
Sect. IV we show and discuss the numerical solutions to the gap and the Bethe- 
Salpeter equations, calculate the Kallen-Lehmann spectral function and analyze 
the particle content of the solutions. Finally in Sect. V we summarize and draw 
conclusions. 

2 The Gaussian effective potential/action 

We shall use the notation and conventions of Ref. [3]. In Sect. II and III we shall 
work in the chiral limit (e = 0), so as to avoid unnecessary complications. Extension 
to the non-chiral case is straightforward. Of course, we use the non-chiral equations 
in the numerical solutions in Sect. IV. 

The effective potential, and the effective action methods in quantum field the- 
ory (QFT) were popularized in the mid-70's. These two are objects with certain 
intriguing theoretical properties: the former is the generating function for the zero- 
momentum one-particle irreducible (OPI) graphs, the latter is the generating func- 
tional for arbitrary momentum OPI graphs [2, 12]. At first only the one-loop pertur- 
bative approximation was calculated in the (j) 4 theory. These two objects need not 
be perturbative, however, and first attempts at their non-perturbative evaluation 
were made slightly later. 

The Gaussian effective potential is a natural product of a variational calculation 
based on the Gaussian ground state trial wave functional [1, 2] of a scalar single 
component cj) quantum field theory (QFT). The O(N) symmetric effective potential 
was calculated e.g. in Ref. [5]. In another original approach Cornwall, Jackiw 
and Tomboulis [CJT], Ref. [11] used certain disconnected ("vacuum") two-particle 
irreducible diagrams to define and calculate the ground state ("vacuum") energy, 
as per Goldstone's (many-body) theorem [13]. The resulting vacuum energy defines 
a (non-perturbative) effective potential [12], which together with the kinetic energy 
defines the effective action. When one minimizes the CJT vacuum energy using 
a particular variational Ansatz, the resulting minimization conditions, or the gap 
equations are equivalent to the ones obtained in the canonical Gaussian variational 
approximation. It is less obvious that the two-body (Bethe-Salpeter) equation in 
the latter formalism is equivalent to the mass equation in the former, in the case of 
spontaneous symmetry breaking. We shall explicate here the proof given in Refs. 
[14,6]. 

Instead of the Cornwall-Jackiw- Tomboulis [CJT] appproach we follow Stevenson, 
Alles and Tarrach's more direct calculation [5] based on the Gaussian vacuum wave 
functionals. The latter authors found the "vacuum" (ground state) energy density 
S(M,n;((f>)) given by Eq. (3.2) in Ref. [3]. By definition [2, 12], the effective 
potential V e s (rrii, ((pi)) is 

VG(mi,{<f>i)) = £{mi, (<&)) - £(rm, (&) = 0). (1) 

One may identify the KI\{mi) term in Eq. (3.2) in Ref. [3] with the familiar "zero- 
point" energy density of a free spinless field of mass mj 1 . This seems to imply that 

Indeed, many finite-temperature/density studies [4] have been based on this observation, as the 
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such one meson states are present in the Gaussian approximation to this theory and 
that ?7ij are their physical masses. That is not the case, however, as we shew in 
Ref. [14, 15] (see, also later): The rrii are merely auxilliary quantities (variational 
parameters) that determine the position of certain particle production thresholds 
and the the corresponding branch cuts, but there are no poles in the propagators 
at those mass values. 

On the other hand it has been shown in two different ways [14, 6, 3] that there are 
three massless (Nambu-Goldstone) and one massive (a) state (with a mass different 
from any of the rrii) in the Gaussian approximation. Thus the physical content of 
the Gaussian approximation, though formally well established, remains one of its 
intuitively most confusing aspects. By working out the connection between various 
formalisms we shall shed more light on this issue. 



2.1 The Gaussian effective potential 

We shall use the fact that the effective potential is the generating function of OPI 
zero-momentum Feynman diagrams. In other words, the n-th derivative of the 
effective potential is the OPI Green's function evaluated at zero external momenta. 
Thus the curvature (the second derivative) of the effective potential with respect to 
the corresponding fields evaluated at the minimum yield the inverse of the two-point 
Green functions T^ 1 evaluated at zero momentum squared p 2 = 0, 

- 



which, in turn, are interpreted as squared masses of the corresponding states by 
way of 

-r« 1 (p = 0) = ml (3) 

thus equating the meson masses with Eqs. (2). It had been tacitly assumed that 
such a definition of a, ir masses 



m 2 



V <^o> 2 J ■ 

\ " 1 / mm 

is equivalent to the "single-particle" , or "gap equation" masses M, fx, i.e. that 
m a = M and m^ = fj,, but that assumption leads to an apparent violation of the 
Nambu-Goldstone [NG] theorem, since n > 0, even in the chiral limit (e = 0). This 
conclusion is incorrect, albeit very common in the literature. Many studies were 
devoted to attempts at its rectification: for example it was shown that \x — ► in the 
1/N — > limit, and that was supposed to restore the NG theorem, see Ref. [5]. A 



zero-point energy and the effective potential point of view offer an "obvious" extension of the zero- 
temperature/density formalism. 
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straightforward evaluation of the derivatives in Eq. (4) yields [6], however, 



2 2 / l + 2Ap [3/mm(0) - yO)] - 24A§J M m(0)/ w (0) ) / Af 2 

^1 - A [3W(0) + 5I MM (0)] + 12A2W(0)/ W (0) J ^ 1 ' 

ml = O^fi 2 , (6) 
where I^(p 2 ), Imm(p 2 ) are logarithmically multivalued functions of p 2 defined by 

M . /■ 1 

mm(Sj V (27r)4[fc2-m2 + k]p-p) 2 -m 2 + ie] ' i j 

where m = M, or fi, explicitly evaluated in Eqs. (3.21), (3.22) and (3.23) in Ref. 
[3]. 

The result Eq. (6), of course, restores the NG theorem, but it does so almost 
per fiat: it gives one no insight into the mechanism that brought it about, e.g. it 
tells us nothing about the sum of Feynman diagrams that leads to it. After all, 
each effective potential generates a certain class of loop diagrams (at zero external 
momenta) that one may wish to identify. It is very difficult to see which class of 
diagrams corresponds to Eqs. (5), (6), and which Schwinger-Dyson [SD] equations 
sum up that class. The first, brief answer to that question was given in Ref. [15]; 
now we shall elaborate on it. Once we have identified the Feynman diagrams one 
immediately sees that Eq. (5) is not quite right. That fact was also recognized 
earlier, Ref. [6], but without reference to Feynman diagrams, and the correction 
was also given there, but only with a formal mathematical justification. Here we 
shall give an explicit Feynman diagrammatic interpretation of that formal definition. 

2.2 The Gaussian effective action 

We shall use the fact that the effective action is the generating functional of one- 
particle irreducible (OPI) Green functions [2], i.e., 

T-Hxy) = ( 62Sm] ) (8) 

where 5 is the functional derivative. Then the correct definition of the particle mass 
is the position of the pole in the two-point Green function, i.e. 



where 




(9) 



and S is the Fourier transform of the (Gaussian) effective action, 

S G m)} = Jd A x(T-V G ), (11) 

T = ^(<9^(</>i)) 2 is the kinetic energy density. Eq. (9) shows the distinction between 
the effective action and the effective potential methods: the effective potential is 
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momentum independent and thus cannot correctly describe the pole in the propa- 
gator, except when the pole happens to be at zero momentum/mass. Thus the mass 
obtained from the effective potential method agrees with the one obtained from the 
effective action only when the mass vanishes. Therefore we must solve only Eq. (9) 
for the mass, as the latter appears on both sides of the equation. This distinction is 
insignificant for massless (NG) states and its importance increases with the mass: 
once the mass crosses the lowest (particle pair) production threshold, it acquires 
an imaginary part that cannot be neglected. Thus the pion mass calculated via 
the effective potential might be OK, because it lies below all hadronic production 
thresholds, but the scalar meson mass is definitely not OK. 

Equations resulting from Eq. (10) in the linear a model have been written down, 
but not solved (except in the trivial NG pion case) in Ref. [6]: Eqs. (10) yield (N 
- 1) = 3 massless states (pions in the chiral limit) and one massive state [a meson) 
whose mass p 2 = m 2 is determined by the roots of the following equation 

_2 M 2(, | ^o[3Imm(p 2 ) + I^(p 2 )-12\Pmm(p 2 )Iw(p 2 )] \ , . 
1 ^ 1 - A Wmm{p 2 ) + 5/ w Cp 2 )] + 12X 2 Imm(p 2 )I^(p 2 ) J 

After making the replacement p 2 = m 2 , and a slight rearrangement Eq. (12) turns 
into 

m 2 M 2 ( 1 + 2A [Mmm {m 2 a ) - I^jml)} - 24X 2 I M M(m 2 a )I l 



(7 



A [3/mmK) + 5/ MM (m2)] + 12X*I MM (ml)I, 




To show equivalence of these results to those of the canonical Gaussian approxima- 
tion, we must first remember that in Ref. [3] we showed the scalar (a) sector coupled 
Bethe-Salpeter (BS) equations that sum infinite classes of connected, though not 
necessarily OPI Feynman diagrams. The mass of the scalar meson is determined by 
the "pole condition" in the scalar channel BS solution, Eq. (3.26) in Ref. [3], that 
reads 



(s-M 2 )V(s) = 0. (14) 



where 

V{s) = 1-A 



M 2 \ I M 2 

3 ( 1 + 3 -^2 )Imm{s)+ 5 + 3— — 2 )l^{s) 



( M 2 \ 

+ 2X 2 I MM (s)I^( S )h+3—j^\, (15) 

is the discriminant of the coupled BS equations (see Ref. [3]). Collecting terms we 
find 



s = M 1 



' \ + 2A [3I M m(s) - I^(s)} - 2AX 2 I MM (s)I^{s) \ 
k I-X [3Imm(s) + 5I^(s)] + 12X 2 Imm(s)I^(s)) 

Upon replacing s = m 2 , this equation becomes identical to Eq. (13) for the a 
mass. In other words, the results of the Gaussian effective action approach are 
exactly identical to those of the Gaussian BS equation (or, in the many-body theory 
language, to the mean-field theory + RPA), proving which was one of our goals. 
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As stated above, in the Gaussian approximation, it is unclear which (classes of) 
OPI diagrams are generated by the effective action. We have readily calculated the 
analytic form of the effective potential curvature at the minimum, Eq. (13), but 
it would be a major challenge to identify the corresponding class of OPI diagrams 
without the benefit of the above insights obtained from the Gaussian Bethe-Salpeter 
equation in Ref. [3]. Furthermore, the factor (s — M 2 ) in the equation (14) explicitly 
shows that the M-particle pole has been "amputated" from the amplitude, i.e., that 
the summed diagrams are one- M-particle irreducible. Some confusion arose due to 
the fact that Tij(p 2 = s) is a two-point OPI Green function, whereas the BS equation 
defines a connected four-point Green function. There is no contradiction, however, 
as one can see after "amputating" the external "legs" of the BS amplitude: the 
result is just an s-channel propagator, i.e., a two-point Green function. We have 
thus given an explicit proof of a formal property of the effective action, but this 
does not begin to tell us what branch of Eq. (13) to solve. 

As the solutions to Eq. (13) are expected to lie in the complex p 2 plane, one must 
specify the sheet ("branch") of the Riemann surface that one is working in. That is, 
in the effective action/potential approach at least, a priori impossible: there is no 
reason why one branch should be preferred to another. Moreover, equations, similar 
to Eq. (13), in models with fermions [16] have been found to contain roots on the 
real axis of the second-, as well as of infinitely many other lower Riemann sheets. In 
the present case there are bound to be even more roots as there are two thresholds 
and two sets of infinitely many Riemann sheets. This ("sheet") ambiguity in the 
effective potential formalism can be resolved by referring to its connection to the 
BS equation. So we turn to the study of analytic properties of the Gaussian BS 
equation. 



3 Analyticity of the Gaussian BS equation 

We shall show the exact equivalence of the solutions to the Gaussian BS equation 
and the so-called N/D equations in S-matrix theory. N/D equations are one way 
of implementing the (two-body) unitarity and causality conditions in a relativistic 
setting, which, in turn translate into analytic properties of the scattering amplitude. 
These analytic properties are important as they tell us what branch of the equation 
to solve and the solution's physical interpretation (bound state, resonance, "anti- 
bound state") [7, 8, 9, 10]. 

The broken-symmetry connected four-point Green function, Eq. (3.12) in Ref. 
[3] , for the scattering of two non-identical ( "pion-sigma" scattering) scalar particles 
has the same generic form of a geometric series as in the symmetric phase, see Ref. 
[1], but with an additional pole term in the "potential" due to the "elementary" 
(massive) "pion" excitation, see Eq. (3.11) and Fig. 4 in Ref. [3]. Such a pole term 
is known in the S-matrix literature as the Castillejo-Dalitz-Dyson (CDD) pole. It 
has been known for some time [17] that such a geometric progression of Feynman 
diagrams corresponds to the solution of the (S-wave) N/D equations in the s-channel 

n / x _ N n(s) ( s 

^"i + ^/^^Wimn^) (i7j 

The solutions to the N/D equations are not unique, however, the arbitrariness show- 
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ing up in the form of so called Castillejo-Dalitz-Dyson (CDD) poles. The position of 
a CDD pole, and the coefficient multiplying it are arbitrary in the usual S-matrix, 
or "bootstrap" approach, but in our approach they are completely determined by 
the Gaussian approximation to the underlying a model Lagrangian. The physical 
interpretation of CDD poles used to be controversial, but the present-day consensus 
is that they correspond to elementary particles/fields in the theory, which conjecture 
is confirmed by our results in the Gaussian approximation. 

We may rewrite the pion-channel kernel ("polarization function") H n (s) of the 
Bethe-Salpeter equation 



^(s) = I Mi i{s) 

in the "dispersive" form (see Eqs. (20) 

d 4 k 1 



(18) 



Imh( s ) 



J (27 



(2vr) 4 [k 2 — M 2 + ie] [(k - P) 2 -n 2 + ie] 



= W°) 



1 

2X t 



(4vr) 2 

2 



Kmll(s) = — 







M 2 




M 2 J (4tt) 2 
Im K Mfl (t) , 



K M u,{s) 



(19) 



where s = P 2 and the real and imaginary parts are 



Im K Mf i(s) 



-Im Imh{s) 



(M - n) 2 



Re K Mli (s) 




(M + n) 2 



6{s-{M + n) 2 ) 



i M 1 L 

log — + - 1 + 
fi 2 \ 



(M - fi) 2 




( M + ^Vanh 



The (momentum) s dependent part of this integral is an analytic function in the cut 
complex s plane. There are in general two logarithmic branch cuts (one stretching 
from M+fj, to +oo, another from M—fi to — oo, though on the first ("physical") sheet 
only the right-hand-side cut appears) determining a Riemann surface with infinitely 
many sheets. This is rather different from the corresponding nonrelativistic case 
which has only one (square root) cut with two sheets. 

Comparing Eq. (3.12) in Ref. [3] with Eq. (17) above, it becomes clear that 
the form of the Gaussian approximation ir propagator demands that the numerator 
function A^-(s) equal the pion channel "potential" V n (s), Eq. (3.11) in Ref. [3] and 
the form of the denominator in Eq. (17): 



N n (s) 



V n (s) = 2A 1 + 



M 2 
i — fi 2 



dt N n (t) 



IT 



t 



If 



Im/ M(U (i) = N n {s) 



W°) + 



7T 



dt 



t{t 



ie 



■Im I M ^(t) 



(21) 



(22) 
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Equation (21) tells us that s = /i 2 is the position, and 2\qM 2 is the coupling 
strength of the CDD pole, whereas Eq. (22) dictates the value of the "subtraction 
constant" 



n*(o) = W°) = 



I (M)-I (fi) 
M 2 — ji 2 




M 2 — e/v \ 
M 2 -fi 2 J 



(23) 



that is, in turn, fixed by the gap equations. 

Since the N/D approximation is unitary by construction we conclude that the s 
channel Gaussian BS scattering amplitude is also unitary. By the same token, one 
can show that the a channel propagator can be written as a solution to the matrix 
N/D equations, see Ref. [8]. The subtraction constants are fixed as in Ref. [3]. 

Now that we have established analytic properties of the Gaussian BS equation, 
we may look for its solutions. As noted above, physically interesting solutions are 
to be found as follows: (1) bound states on the real axis of the physical sheet, below 
all thresholds; (2) resonances in the fourth quadrant of the "second" sheet; (3) "an- 
tibound states" on the real axis of the second sheet, below both thresholds. There 
are also new kinds of S-matrix singularities in the relativistic quantum field theory 
(QFT) that do not appear in nonrelativistic quantum mechanics. One such possible 
new relativistic singularity is the CDD "pole" (or "zero" in the denominator). CDD 
poles are associated with "elementary" particles/fields in the theory, see p. 400 in 
Ref. [7]. In our case this means a field in the a model Lagrangian, If we remember 
that the linear a model has been shown to be the low-energy limit of the Nambu- 
Jona-Lasinio (NJL) chiral quark model, we may say that the "elementary" particles 
here are just the bare NJL quark-antiquark states. These states lead to real poles 
only in the weak meson-meson coupling limit, as we shall show below. Thus the a 
resonance may become a stable state only in the limit of weak interactions. 

It should be noted, however, that Imm(s), I^n(s) are analytic functions with 
imaginary parts above the corresponding thresholds, so that Eq. (14) actually 
implies two equations: one for the real and one for the imaginary part. Usually 
only the real part is considered, however. That is all right if the (real part of the) 
root lies below all thresholds. If it does not, as in our case, one must look at the 
equation for the imaginary part, as well. As one moves away from the real s axis, 
each equation yields a line of roots in the complex s plane. The intersection of the 
two lines (the real and imaginary roots) then yields the position of the pole. Only 
the pole on the second lower sheet 2 , if it exists at all, determines the mass and 
width of the resonance. 

We have numerically solved the real part of the scalar meson mass/Gaussian 
BS equation on the real axis of the physical sheet, looking for bound states. The 
results are shown in Fig. 9 of Ref. [3] . Note the double- valuedness of the solutions. 
The question arises: can one ascribe particle-like properties to the heavy branch of 
the solution. Short of a pole search in the second sheet, that question can only be 
answered by calculating the Kallen-Lehmann spectral function. 

2 Imm{s), Iu,u(s) have logarithmic branch points and therefore infinitely many sheets, in contrast with 
the nonrelativistic case where the branch points are of the square root type, with only two sheets. 
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4 Kallen-Lehmann Spectral function 



As argued in Sect. Ill, the Gaussian pion propagator D n (s) is an analytic function 
in the cut s-plane and, as such, it allows a dispersive, or Kallen-Lehmann represen- 
tation 

ArW = " / dt - pAt) , , (24) 
J t — s — te 

where 

Ptt(s) = --lmD^(s) , 

7T 

is the spectral density function. The latter represents the mass distribution of 
physical excitations in this channel. In Ref. [15] we have explicitly shown that in 
the pion (J p = _ ) channel and in the chiral limit, the spectral function 

Pn (s) = aS(s) + c(s)9(s — (M + fi) 2 ) , (25) 

contains only one Dirac delta function, instead of two, as naively expected. In 
Fig. 1 we show the same spectral density in the nonchiral case (explicitly broken 
0(4) symmetry) and again find only one Dirac delta function, this time at yfs = 
140 MeV. Another way of saying this is that the strength with which the state at 
s = fi 2 appears in the spectrum is zero, i.e., the state decouples from the single- 
particle spectrum. The heavier excitations 4>i,i = 1,2,3 correspond to unstable 
quasi-particles [18] that decay into an odd number of lighter Goldstone bosons. Thus 
there is no particle doubling in the Gaussian approximation, contrary to suggestions 
e.g. by Tornqvist [19]. 

Similar comments hold for the a sector, the analysis being more complicated 
due to two different kinds of intermediate state being possible there. The sigma 
channel is phenomenologically more interesting than the pion one because that is 
where many experimental "supernumerary" states have been observed. In Figs. 2, 
3, and 4 we show the scalar Kallen-Lehmann spectral density 

Pa(s) = --lm(D a (s)) , 

7T 

at various free parameter values. Here the scalar propagator is defined as 

o \ 

= ( S -M>)W <26> 

the remaining terms in the Dij matrix elements contribute to the effective 0"^^, i = 
0,1,2,3 vertex form factors. Once again we find only one bump, or Dirac delta 
function in the density of states, depending on the coupling strength Ao and other 
free parameters. 

We showed in Fig. 8 of Ref. [3] that the scalar channel (dressed a meson) pole 
position is always shifted downward from the "elementary" sigma field's (0°) CDD 
pole at s = M 2 , in accord with the variational nature of the Gaussian approxima- 
tion. As the coupling constant Ao drops below some critical value A c , which is a 
function of the masses fi, M and the cutoff A, the a meson mass becomes purely 
real as its position drops below the 4/j 2 threshold, see Fig. 2. Thus, within the 
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Gaussian approximation, the a meson can be (quasi-)stable at weak couplings, but 
at moderate and strong couplings it is always a broad resonance. The behaviour 
of the critical value of the a meson mass m° c as a function of the cutoff A can be 
gleaned from Fig. 8 in Ref. [3] (and similar graphs for values of A between 0.4 and 
1 GeV): the crossing point of the m cr (M) and 2/j,(M) curves lies around M ~ 300 - 
400 MeV until it disappears altogether for A > 0.5 GeV. 

In Ref. [3] (see Fig. 9.) we showed that at intermediate values of the coupling 
Ao and low cutoff A there is a second, much heavier solution to the real part of the 
BS equation besides the usual light one. Thus there is once again the possibility of a 
second pole in the S matrix. Such a conclusion would be premature, as can be seen in 
Fig. 3, and 4: there is no enhancement of the scalar spectral function at high energy, 
whereas at low s the standard solution pole can be seen as a Dirac delta function 
turning first into a narrow then a wider peak as one moves up in energy. Thus 
we see that there is no enhancement in the density of states at the corresponding 
energy/mass and thus there is no second resonance peak. We conclude that the 
heavy solution exists only in the real part of the BS equation, whereas the imaginary 
part does not have a root in the vicinity. 

5 Summary and Conclusions 

The mean-field method was initially fraught with problems when applied to the 
linear a model with spontaneously broken internal symmetry - the Goldstone the- 
orem did not seem to "work". This problem was solved, at the price of opening 
new questions [14]: The Goldstone boson found in the Gaussian approximation [14] 
turned out to be a composite state that apparently co-exists with the massive "ele- 
mentary" state with identical quantum numbers. A similar situation occurs in the 
scalar sector. At first this looks like a "doubling" of the number of scalar states, 
sometimes invoked in the phenomenological literature on the a meson [19]. We have 
showed here that this doubling is only apparent: we investigated the question of the 
particle content of the (j) 4 theory in the Gaussian approximation by employing the 
Kallen-Lehmann representation [2]. Thus we found that the massive "elementary" 
pion does not appear in the Kallen-Lehmann spectral function as an excitation in 
the pseudoscalar channel. This is the same as saying that the Castillejo-Dalitz- 
Dyson (CDD) "pole" in the ir channel has a vanishing residue in the MFA to the 
linear sigma model. Similarly, only one state has been found to exist in the a 
channel of the MFA. Moreover, we showed that the mass of the composite a state 
agrees exactly with that calculated in the CJT formalism. Another interpretation, 
(that might be only semantically different from the above one) was derived by using 
operator many-body ( "quasi-particle RPA") methods [18, 13]. In that theoretical 
framework the "lighter" states are particles and the "heavier" ones are so-called 
quasi-particles that appear due to the interactions of the particles. 

In summary, we have: 1) established equivalence of the standard Gaussian ap- 
proximation with the effective potential method; 2) shown equivalence to the N/D 
equations and thus proven unitarity. 3) solved the resulting equations and investi- 
gated the theories particle content via the Kallen-Lehmann representation. 

One of the authors [V.D.] wishes to thank Prof. Paul Stevenson for valuable 
conversations and correspondence relating to the Gaussian approximation in general 
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and some of the problems treated here in particular. The same author would also 
like to acknowledge a COE Professorship at RCNP for the year 2000/1 at which 
time most of this work was done. 

References 

[1] T. Barnes and G.I. Ghandour, Phys. Rev. D 22, 924 (1980); see also eds. 
L. Polley and D.E.L. Pottinger, Variational Calculations in Quantum Field 
Theory, World Scientific, Singapore (1987). 

[2] R.J. Rivers Path Integral Methods in Quantum Field Theory, Cambridge Uni- 
versity Press, Cambridge (1987), and B. Hatfield, Quantum Field Theory of 
Point Particles and Strings, Addison- Wesley, Reading (1992). 

[3] Issei Nakamura and V. Dmitrasinovic, Prog. Theor. Phys. 106, 1195 (2001). 

[4] G. Amelino-Camelia, J. D. Bjorken, and S. Larsson, Phys. Rev. D 56 6942 

(1997) ; S. Chiku and T. Hatsuda, Phys. Rev. D 57, R 6 (1998), and Phys. 
Rev. D 58, 076001, (1998); H.-S. Roh, and T. Matsui, Eur. Phys. J. A 1, 205 

(1998) ; J. T. Lenaghan, D. H. Rischke, J. Schaffner-Bielich, Phys.Rev.D 62, 
085008 (2000); Y. Nemoto, K. Naito and M. Oka, Eur. Phys. J. A 9, 245 (2000). 

[5] P. M. Stevenson, B. Alles and R. Tarrach, Phys. Rev. D 35, 2407 (1987). 

[6] A. Okopihska, Phys. Lett. B 375, 213 (1996). 

[7] S. Gaziorowicz, Elementary Particle Physics, (J. Wiley and Sons, New York, 
1966). 

[8] K. Nishijima, Fields and Particles, (W.A. Benjamin, Reading, MA, 1969). 

[9] S. C. Frautschi, Regge Poles and S-matrix Theory, (W.A. Benjamin, New York, 
1963). 

[10] P. D. B. Collins and E. J. Squires, Regge Poles in Particle Physics, (Springer, 
New York, 1968). 

[11] J.M. Cornwall, R. Jackiw and E. Tomboulis, Phys. Rev. D 10 2428 (1974). 

[12] C. Itzykson and J.B. Zuber, Quantum Field Theory, (McGraw-Hill, New York 
1980). 

[13] A. Fetter and J.D. Walecka, Quantum Theory of Many-Body Systems, 
(McGraw-Hill, New York 1971). 

[14] V. Dmitrasinovic, J. A. McNeil and J. Shepard, Z. Phys. C 69, 359 (1996). 

[15] V. Dmitrasinovic, Phys. Lett. B 433, 362 (1998). 

[16] V. Dmitrasinovic, Phys. Lett. B 451, 170 (1999). 

[17] F. Zachariasen, p. 97 - 99 in Recent Developments in Particle Physics, editor 
M. J. Moravcsik, (Gordon and Breach, New York, 1966). 

[18] H. W. L. Aouissat, O. Bohr, and J. Wambach, Mod. Phys. Lett. A 13, 1827 
(1998). 

[19] N. Tornqvist, Phys. Lett. B 426, 105 (1998). 



12 



Figure 1: The isovector pseudoscalar (pion channel) spectral density p n (s) in the Gaussian 
approximation to the 0(4) linear sigma model as a function of the CM energy y/s for 
various values of the free parameters. The vertical line at y/s = 140 MeV represents the 
Dirac delta function. Note the absence of other delta functions. The threshold of the 
continuum is at y/s — M + p,. 



Figure 2: The isoscalar scalar (sigma channel) spectral density p a (s) in the Gaussian 
approximation to the 0(4) linear sigma model as a function of the CM energy y/s for A 
= 1 GeV and various values of the variational parameter (mass) M, i.e. of the coupling 
constant Ao- The vertical line represents the Dirac delta function. The first threshold of 
the continuum is at y/s = 2p, the second threshold (2 "a") is at y/s = 2M where a cusp 
in the spectral density can be seen. In both cases the physical a mass lies below the lower 
threshold at y/s = 2/i, so the a meson is stable. 



Figure 3: As in Fig. 2, but for cutoff A= 0.4 GeV and three values of the variational 
parameter M = 330, 370 and 500 MeV (with p— 146, 156, and 165 MeV, respectively). 
In one case the physical a mass lies below the lower threshold at y/s = 2/i, in the other 
two it lies above their respective thresholds, so accordingly the a meson is either stable 
or a resonance. Note the widening of the a resonance peak as its mass increases, and the 
coupling constant grows Aq with it. 



Figure 4: Same as in Fig. 3, but rescaled to show higher values of y/s: note the absence 
of any enhancement in the region of the "second solution" (the cusp is due to the second 
threshold). 
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